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How entangled is a randomly chosen bipartite stabilizer state? We show that if the number of 
qubits each party holds is large the state will be close to maximally entangled with probability 
exponentially close to one. We provide a similar tight characterization of the entanglement present 
in the maximally mixed state of a randomly chosen stabilizer code. Finally, we show that typically 
very few GHZ states can be extracted from a random multipartite stabilizer state via local unitary 
operations. Our main tool is a new concentration inequality which bounds deviations from the mean 
of random variables which are naturally defined on the Clifford group. 

PACS numbers: 03.67.Mn, 03.67.-a 

Keywords: entanglement, stabilizer formalism, concentration of measure 



I. INTRODUCTION 

Randomly chosen states and subspaces play a central 
role in the study of quantum information. For exam- 
ple,the consideration of random stabilizer codes played 
a crucial role in one of the first proofs that there exist 
good quantum correcting codes [l5|, as well as much of 
the early understanding of entanglement distillation and 
quantum channel capacities [3, [Til, [23. More recently, 
through an improved understanding of the typical prop- 
erties of randomly chosen quantum states, ex pre ssions for 
many capacities of quantum channels 3i [9|,127|. [33| . sev- 
eral advances in cryptography [J, [TtI . 18l |. and an emerg- 
ing understanding of quantum correlations in high di- 
mensional systems [l9j have all been attained. 
Perhaps the property of quantum states that it is most 
important to understand is entanglement. Entanglement 
is an essential resource in quantum information which 
nevertheless remains quite poorly understood in general. 
Even in the asymptotic limit it is difficult to character- 
ize the entanglement in a bipartite mixed state. Indeed, 
such fundamental quantities as the entanglement of for- 
mation, Ep, and the distillible entanglement. En, are 
unknown in all but a few examples (see, e.g..|24l.[33.[34|). 
It has, however, proved possible to find tight bounds 
on the typical Ep and Ep) of a random mixed state, 
Pu = ^ U\i){i\W of rank 2^, both local dimensions 

roughly 2" 3> 1, and where U is distributed according 
to the unitarily invariant measure (i.e., the Haar mea- 
sure) on ILJf(i4dR)[l9|. The surprising result is that with 
high probability, pu has Ep « n and Ep) < n — ^ which 
implies that a typical pu of high rank (2*^ « 2^") has 
near maximal entanglement of formation while having 
distillable entanglement which is exponentially smaller, 
implying cither an extreme irreversibility in the creation 
of Pu or a near maximal violation of the conjecture of 
Ep^s additivity [2^. While this dichotomy is quite strik- 
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ing, its physical and computational significance arc not 
at all clear - no constructions of such extreme states are 
known, and the generation of a state distributed like pu 
would require exponential resources. Characterizing the 
typical entanglement of random states whose distribu- 
tion can be generated efficiently is thus crucial to under- 
standing whether this irreversibility is a fact of nature or 
merely a mathematical curiosity. 

In this paper, we characterize the typical entanglement 
for just such a distribution - the uniform distribution on 
the set of stabilizer states, which can be generated effi- 
ciently using the random walk based algorithm of [lo| . 
Stabilizer states are relevant to almost all known quan- 
tum error correcting codes, and as such it is hoped that 
a characterization of their typical entanglement proper- 
ties will not only shed light on the irreversibility question 
mentioned above, but also point us towards better codes. 
Furthermore, since a highly entangled multipartite sta- 
bilizer state is the fundamental resource in the one-way 
model of quantum computation [2^, a deeper under- 
standing of such states may elucidate the role played by 
entanglement in quantum computations. 
The bipartite entanglement of stabilizer states has previ- 
ously been explored in [l], |3, [3l , with the result that any 
such state can be transformed by local unitarics (i.e., it 
is LU-equivalent) to a tensor product of EPR pairs and 
(possibly classically correlated) local states. An expres- 
sion for the number of EPR pairs that can be extracted 
from any particular stablizer state in terms of the struc- 
ture of its stabilizer group was also found. Our contribu- 
tion is to estimate the expectation of this expression for 
a random stabilizer group, and provide an exponential 
bound on deviations from this estimate. Shortly after 
this work first appeared, results for the case of pure sta- 
bilizer states were presented in Q . 

In contrast to a rank 2*^ state with a Haar-induced distri- 
bution, wc find that a random 2" x 2" rank 2*^ stabilizer 
state has Ep = Ep) ^ n — ^. That Epi and Ep must co- 
incide is clear, since any stabilizer state is LU-equivalent 
to the tensor product of EPR pairs and a separable state. 
It is facinating, however, that the value they take is es- 
sentially maximal - a stabilizer state of rank 2^ can have 
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at most 2n — k pure qubits, and in a typical such state 
each of these is half of an EPR pair. 
In addition to these results on bipartite entanglement, 
we are able to characterize the typical number of GHZ 
states that can be extracted via local unitaries from a 
random multipartite stabilizer state. Unlike EPR pairs, 
which are abundant in a random bipartite state, GHZ- 
equivalent states are quite uncommon. For example, we 
find that for a pure m-partite state in which all m sys- 
tems are of size roughly n qubits, in the limit of large n 
the expected number of GHZ states that can be extracted 
via local unitaries is close to zero, and that significant de- 
viations from this mean are unlikely. 
We will concentrate on the case where all of our sys- 
tems have asymptotically equal numbers of qubits (e.g., 
bipartite systems of ua + ns qubits with ha = n and 
ub = n + 0{\ogn)), both because this case is likely the 
most useful in terms of applications and because it is 
exactly where standard Markov-type arguments break 
down. For instance, in the case of bipartite stabilizer 
states a straightforward Markov inequality argument can 
be used to give bounds on deviations of entanglement 



from its mean of the form 



but when ua = n and 



Ub = n + O(logn) this bound is quite weak, scaling like 
^oV) • Using a new concentration inequality, we are able 
to provide exponential bounds in this regime. 
Our main tool throughout is Theorem IH.ll which cap- 
tures the notion of measure concentration on the Clifford 
group. In particular, this Theorem is a quantitative ver- 
sion of the intuitively obvious observation that a slowly 
varying function on the Clifford group won't deviate sig- 
nificantly from its mean. The result is quite general, and 
we expect it will prove useful in further analyses of the 
entanglement of stabilizer states, as well as the analysis 
of stabilizer codes. 

The paper is organized as follows. In Section |TT] we 
present a bound on deviations from the mean of random 
variables on the Clifford group. In Section IIIII we use 
this inequality to characterize the typical entanglement 
in a pure bipartite stabilizer state, in section HVl we study 
pure multipartite states, while in Section |V] we turn our 
attention to mixed bipartite stabilizer states. Section IVll 
contains a few comments on other applications of our in- 
equality as well as some open questions. 
We use the following conventions throughout, log and 
exp are always base 2. The Pauli group on n qubits is 
denoted by Pn- An abelian subgroup of Vn with 2""'^ 
elements will typically be called S"^'' and have gener- 
ators {Si}^~^. Two elements of Vn, Pi and P2, either 
commute or anticommute with the commutation relation 
P1P2 = serving as a definition for uj{Pi,P2). 

Angle brackets will denote the group generated by the el- 
ements they enclose, so that, e.g., S'""'' = {Sj^'l^i- The 
dimension of a subgroup of Vn is the logarithm of the 
number of elements in the group, so that dimS'"^'^ = 
n—k. We say that j?/') is stabilized by U when U\il)) = ip, 
and call \ip) a stabilizer state on n qubits if it is simulata- 
neously stabilized by all elements of a maximal Abelian 



subgroup (S) of the Pauli group on n qubits. A mixed 
stabilizer state of rank 2'' is the maximally mixed state on 
the subspace stabilized by an Abelian subgroup S C Vn 
of size 2""*^, or equivalently is the maximally mixed state 
on the stabilizer code defined by S. 

The Clifford group on n qubits is denoted by C„ and 
its elements arc typically called c. A real valued func- 
tion, F, on a metric space {X, d) is called 77-Lipschitz if 
\F{x) — F{y)\ < r}d{x,y) for all x,y G X. E denotes an 
expectation value, while g C^r G is a random variable 
distributed uniformly on G. 



II. A CONCENTRATION INEQUALITY ON 
THE CLIFFORD GROUP 



The notion of measure concentration is a generalization 
of the basic fact from probability theory that the empir- 
ical mean of many i.i.d. random variables, j^'^iLi -^i, 
tends to be very close to the mean of the underlying 
distribution. The point is that not only the empirical 
mean of a large number of random variables, but any 
function which depends in a sufficiently smooth way on 
a large number of fairly independent random variables 
will tend to be roughly constant. The imprecision of of 
the previous sentence is a reflection of the broad range of 
problems this idea can be applied to - one's definition of 
"fairly independent" or "roughly constant" dep ends on 
the particular question under consideration [20l. [sol. [Slj . 
To make precise the notion of a smooth function on the 
Clifford group, we must first introduce a notion of dis- 
tance between two elements of the group. One natural 
candidate is an analogue of the Hamming distance on the 
set of binary strings. In particular, for some fixed set of 
generators of Vn, {Si}-fl^i, we let the {S'i}-distance be- 
tween ci , C2 e Cn be the number of generators on which 
ci and C2 disagree (ignoring differences in phase). 



}(C1,C2) = #<i 



and choose the smallest such value over all generating 
sets of Vn'- 



d{ci,C2) = mind{5^}(ci,C2). 

{Si} 



(1) 



That this defines a metric is shown in section HVl 
Our "smooth" functions will be those which are 1- 
Lipschitz. That is, we will study deviations from the 
mean of real functions F on C„ such that \F{ci) — 
F{c2)\ < d(ci,C2). The precise meaning of the claim 
that they are "roughly constant" is given by the follow- 
ing theorem. 



Theorem II. 1 Let 

{Cn, d) and c C„ 
variable. Then 



F be a 1-Lipschitz function on 
be a uniformly distributed random 



P{\F{c) ~EF{c)\ >5) < 2 exp 



64n 
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We will prove this theorem, which is quite similar to a 
result of Maurey for the symmetric group [l^l , by using a 
result of [13] that characterizes concentration on a finite 
metric space. In particular, we say that a metric space 
{X, dx) has length at most L if there exists an increasing 
sequence of partitions (i.e, a filtration) of X , 



{X} = x° c • • • C x" = 



ex 



and real numbers ao,...,a,„ with J^i^Qof — LP' such 
that if = {^j}j=i...r. and A^^,A\. C ^j,"^ there exists a 



bijection (/i^^ : —* A\. that satisfies dx{x, (l>lp {x)) < ai 
(Fig. ^ may help elucidate this fairly clumsy definition). 
A concentration inequality is then given by the following. 



Theorem II.2 (^.2 of fM]) Let {X,dx) be a finite met- 
ric space of length at most L and let P be the normalized 
counting measure on X. Then, for every 1-Lipschitz func- 
tion F on {X, dx) and every S > 

P {{F >EF + S}) < exp [-'S^/2L^] . 



Theorem III. II is an immediate consequence of the follow- 
ing theorem. 

Theorem II. 3 The length of {Cn,d) is at most \/32n. 

Proof Let {Si}^2i be a set of generators for Vn such 
that for t G {1, . . . , n}, S2t-i and 5*24 anticommute and 
all other pairs of generators commute. We choose our 
filtration of Cn to be 

where {P^j^^i are independent elements of Vn such that 
ujiP^,Pj)=u;{S„Sj) and 



P1P2 



{c e Cn\cS^c^ ^ P,hr i=l...k}. 



Now we need to find bijections 



^PlP2...Pk-lP 



^PlP2...Pk-lQ- 



Wc first consider the case where Q is independent of 
{Pi, . . . , Pfc_i, P} and k is even. Since both P and Q 
anticommute with Pk-i and commute with {Pi\iZi^ '^e 
can always find Ti and T2 such that the ordered lists 



and 



{Pi,...,Pfc_i,p,Q,ri} 



{Pi,...,Pfc„i,g,p,r2} 



have the same commutation relations and are indepen- 
dent generators for the same group. In particular, given 
any Ti which commutes with {Pi\^Zi and is independent 
of 

{Pi,...,Pfe_i,P,0} 



we can simply choose 

r2 = Ti 

\iuj{P,Ti) = w(Q,Ti) and 

= Pfe-iPi 

if ^(P, Pi) 7^ Pi). This allows us to extend both sets 
to generators of P„ with the same Pk-^^s, ■ ■ ■ P2n - That is, 

~ {Sfl-i = {Pi, . . . , Pk-l,P, Q, Pi, Pfc+3, • ■ • , P2n} 

and 

S'^ = {S?h = { A, . . . , Pk-l, Q, P, P2, Pk+3, . . . , P2n} 

both generate the Pauli group and have the same com- 
mutation relations. As a result, there is an element of 
the Clifford group cpg e C„ such that 



and choosing 



CPQ^i CpQ = 



^Pl''..P,._l(ci) = ClCpQ 



gives us 



d{ci,(t)pf p^_^{ci)) < f^{ctsfci}(ci>0Pi'^..p,._i(ci)) 
If k is odd and Q is independent of 
{Pi,...,Pfc_i,P} 



PQ 



< 3. 



we can make a similar argument. In particular, if PQ 
—QP (i.e., a;(P, Q) = 1), we can immediately extend 



and 



{Pi,...,Pfc_i,p,g} 
{Pi,...,Pfc_i,0,p} 



to 



•5"^ = {Sih = {A, . . . , Pfc-l, P, 0, Pfc+2, . . . , P2„} 



and 



S'^ = {S^h ^{Pl,---, Pk-l, Q, P, Pfc+2, . . . , P2„}, 

so that and S'^ have the same commutation relations 
and are different in only two entries. Choosing 



where again 



''p?...Pk-Mi) = ci(^PQ^ 
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wc find 

If P and Q commute, tliere are Ti and T2 sucli that 
{Pu...,Pk-i.P,Q,Ti,T2} 
are independent, and satisfy 







QTi 




QT2 


= -TiQ. 


PTi 


= -TiP, 


PT2 


= T2P, 



so that there are 

S"" = {Sfh = {Pi, . . . , Pfc^l, P, Q, Ti, T2, Pfc+4 . . . , P2„} 

and 

= {5?}. = {Pl,...,Pfc-l,0,P,T2,ri,Pfc+4,...,P2n} 

with the same commutation relations. Once again we use 
cpQ such that 

cpqS[c^pq = Sf 
to define tpp^ p^__^ and find that this time, 

rf(ci,0P^^..P,_,(ci)) <4. 

When Q e (Pi, . . . , Pfe_i, P) and fc is even, the re- 
quirements that Q ^ (Pi, . . . , Pfc_i), ct;(Pi,Q)=0 for 
i=l,...,fc— 2, and uj{Pk-i-, Q)=l imply that the only 
choice for Q that is not equal to P is just Q cx PPt-i, 
so we can let 



— {Pi , • ■ • , Pfc-l , P, Pfc+l • ■ • , P2n} 



and 



SQ = {S^}, = {Pi, . . . , Pfe_i, Q, Pfc+i, . . . , P2„} 

and proceed as above, with the result that 

rf(ci,0P^^..P,_,(ci))<l. 

If k is odd and Q G (Pi, . . . , Pk-i,P), we find that Q 
can only satisfy the required commutation relations and 
belong to (Pi , . . . , Pfc-i , P) if Q cx P, so in this case we 
find that 








X 
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FIG. 1: Finding the length of Cn- We construct an increasing 
sequence of partitions of C„, {Cn} = ^ ' ' ' ^ X^" = 
{{c}}cgc„ , such that any pair of sets in a partition which 
belong to the same set in the preceding partition have a bi- 
jection, <j}, that satisfies d{c, (jy^c)) < 4. Our filtration contains 
2n partitions, so our length is then V32n. 



(recalling that we equate Clifford group elements which 
differ only by phases). 

Collecting the various cases, we see that we can find the 
bijections we require with ak < 4, so that the length of 
Cn is no more than 



L = 



\ 



2n 



^ al =. V32n. 

fe=i 



□ 



III. PURE BIPARTITE STATES 

Recently, studied the entanglement of a bipartite 
stabilizer state in terms of the structure of its stabilizer 
group. Their result can be summarized as follows. 

Theorem III.l (Result 1 of fZ3 /) Let IV'ab) be a pure 
bipartite stabilizer state with stabilizer S. Then, [ip) is 
LU-equivalent to 

E{\iI^ab)) = ^ [dimS'-dim(S'^+5^)] = n^-dimS'^ 

(2) 

{.9 6 S\g = gA(S) Ib}- 



EPR pairs, where = {g G S\g ~ Ia® gu} and Sg = 



Using this theorem, we will investigate the average en- 
tanglement of a stabilizer state, then strengthen our re- 
sults to statements about typical states using Theorem 
III. II We begin with the following lower bound. 
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Theorem III. 2 Let \'ip) he uniformly distributed on the 
set of stabilizer states on AB , where A contains qubits 
and B has ub qubits and ua > ns- Then 



nsiiJA)] > UB - 



2nA ' 



where ipA = T^'Cb \^Jab){'>1jab\ and S'(V'a) = 
— Tr ipA log ipA is the von Neumann entropy of ipA ■ 

Proof We use a result of [2l| for the average subsystem 
purity of a state uniformly distributed over the entire 
Hilbcrt space AB, together with a result of which 
implies that the average over stabilizer states takes the 
same value. In particular, in [ssj it was shown that if \ip) 
is a uniformly distributed pure state on AB, 



2nA+nB + 1' 



where the expectation is with respect to the uniform mea- 
sure on all states in AB. Furthermore, the observation of 
(lol | that a so-called bilateral Clifford twirl is equivalent 
to a bilateral full twirl implies that the average purity of 
a random stabilizer state has the same value. That is. 



2nA+nB -\- 1 



where the expectation is with respect to the uniform dis- 
tribution on stabilizer states. To complete the proof, we 
use the fact that — logTrp^ < S{p) together with the 
concavity of the log function to conclude 



ESiM > -ElogTrV'^ 
> -logETrV'A 



> log 



> Ub 



2"-4 - 

2"-4 ■ 



2"-E 



□ 

We will also need the following lemma, whose proof 
depends on a more general lemma of Section IIVI 



Lemma III. 3 As a function of c £ CnA+nsj entan- 
glement of \^) = cjO)'*''"'*"'""^-' is 1-Lipschitz with respect 
to the metric defined in Eq (QJ). 



Proof Using Eq. ([2|) together with Lcmma riV.4[ which 
is proved below, immediately implies the result. □ 
Theorem IIIL2| which estimates the average entangle- 
ment of a bipartite stabilizer state, can be combined with 
this evaluation of the Lipschitz constant of the bipartite 
entanglement of stabilizer state li/)) = c|0)'^("^+"«) to 
yield a characterization of the typical entanglement in 
such a state. That is, we can use these to prove Theorem 

MM 



Theorem III. 4 Let \tjj) be uniformly distributed on the 
set of stabilizer states on AB , where A contains ua qubits 
and B has ub qubits and ha > ub. Then the probability 
of the entanglement of \ip) deviating from its mean is 
given by 



P iSiiJA) < nSiM] -S)< cxp 



64 (ha + ub) 



where '4>a = Tr^ I-^abX^/'ab | and S{iPa) = 
— Tr ipA log ipA is the von Neumann entropy of tpA ■ 
In particular, letting ua = n + alogn > n = ub, S = ne, 
and considering n >2/e leads to 



P(5(Va) <n(l-e)) <cxp 



2n 



512 2n + alogn 

Proof of Theorem IIII.4I From Theorem III. II and 
Lemma IlII. 31 we can immediately conclude that 

5^ 



P{S{4'a)<^S{iPa)- 5) <c^p 



6A{nA + Ub) 



From Theorem llIL2] wc know that E[S{iPa)] >nB-^, 
so that 

P ( S(tpA) <n - — - ne/2 I < exp 
V n" J 



256 2n -|- alogn 



which leads to 

P{S{^a) <n{l~e)) <cx_p 



ne 



2n 



512 2n + alogn 



□ 



IV. PURE MULTIPARTITE STATES 

The results of the previous section have immediate con- 
sequences for the number of GHZ states that can be LU- 
extracted from a random stabilizer state. In particular, 
we find the following theorem. 

Theorem IV. 1 Let iVAn) be a state uniformly dis- 
tributed on the set of pure m-partite stabilizer states with 
each party holding n qubits, and where m > A. Then if 
for every A < m! < m we let A™ (|'0m)) denote the maxi- 
mal number of m' -GHZ states, ^(|0)®™' + |1)«^™'), that 

can be extracted from IV'm) via unitaries which act lo- 
cally ( with respect to a partition of the m parties into m' 
groups ), 

p(a'"(|V;™)) >enl 



< 



exp 



-[m/2\n—il-l/[m/2]f{l-l/m) 



and 



P{ A'" (IVAn)) >£« <CXp 



64 TO 
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Proof We first consider the case where m! = m. Let 
k = [m/2j, B denote parties 1 through k and A de- 
note the rest, and consider the number of EPR pairs 
with respect to the A\B partition that can be extracted 
from IV'm) via local unitaries on A and B. Supposing 
A™ (I ■0m)) =5, we can see that the total number of EPR 
pairs that can be LU-extractcd between A and B is no 
larger than kn — {k~ l)g as follows. First notice that lo- 
cal unitaries do not alter the local entropies of A and B, 
so that the number of EPR pairs A and B can extract 
is no more than S{B). However, the kg qubit support 
of B's part of the g GHZ states contains only g bits of 
entropy, since the reduced state on these qubits is of the 
form (i|0)(Op'^ + i|l)(l|®'^)®^. B's remaining kn - kg 
qubits can have a maximum of fcn — kg bits of entropy, 
leading us to conclude that S{B) < kn — kg + g. 

Letting a — {m — 2[m/2j )(n/logn) and using Theo- 
rem |nL4] we thus find the probability that A™(|-0„i)) > 
en is no larger than 



P 



{Si^ps) < fcn-(fc-l)ne) 



< exp 



kne^ , , ^„2kn 
512 ^ ' ' mn 



Similarly, for 4 < m' < m, let B denote the three 
smallest groups in our partition and let k be the number 
of parties in B. Once again, if A™ (iV'm)) > the 
number of EPR pairs that can be extracted between A 
and B can be no larger than kn—{k—l)ne. Thcorcm lHL4l 
can then be used, this time with a = (m — 2k){n/ logn), 
to show that 



P {S{ipB)<kn-{k-l)ne) < exp 



fcrte^ , , , 2kn 

-—-{1-1/ky — 

512 mn 



The expressions in the theorem are obtained by sub- 
stituting the values for k and in the case of m' < m using 
the fact that fc > 3. □ 

We cannot understand the entanglement of a tripar- 
tite stabilizer state by simply considering the bipartite 
entanglement of various partitions. In this case, we use 
following theorem, which was proved in Q. 

Theorem IV.2 (Theorem 3, Corollary 2 of Q/j Let |V') 
be a pure m-partite stabilizer state with stabilizer S. 
Then the number of GHZ states, |1'+) = ^(|0)^™ -I- 
jl)*^™), extractable from Itjj) via local unitaries is 

AiS) = dim{S)-dim{Sioc), (3) 

where Sioc is given by Sioc = J2'a=i and Sa ^ {g & 
S\g acts trivially on a}. 

Below we will find an upper bound for the expected 
value of Eq. ([3]) when m — 3, which wc will then combine 
with the following lemmas, which imply that the number 
of GHZ states LU-extractable from a random m-partite 
stabilizer state (to fixed) concentrates tightly around its 
mean value when the number of qubits each party holds 
is large. 



Lemma IV. 3 Consider the binary representation ofVn 
on F2", wherein c G Qn is represented by an element of 
GL2n{¥2) (see, e.g., ]A])- The metric of Eq (0]) is given 
by 

d{ci, C2) = 2n — dimKer[c'i — 62], 

where, e.g., c £ GL2,i(F2) is the representative of c € Cn. 

This lemma, together with the fact that Rank(c'i —63) < 
Rank(c'i —62) + Rank(c'2 — C3), makes clear that the dis- 
tance defined in Eq. [l]is in fact a metric. 
Proof (of Lemma IV. 3) To see this, note that 



min#{i|ciS'ic{ ^ C2Sicl} 

{Si} 

min#{i|(ci - 62)81 ^ 0} 

{Si} 

2n — dim Ker[ci — 62]. 



(4) 



□ 



Lemma IV. 4 As a function 0/ c G C„, where n = 
X]q=i "^c; dimension of Sf^^ (defined in Theorem 
\IV.S^) of S'^ = cSqc^ for some fixed stabilizer Sq is m- 
Lipschitz with respect to the metric defined in Eq (Qp. 
In particular, the number of m-partite CHZ states (with 
m>3) that can be LU- extracted from the state with sta- 
bilizer S'^ = cSoc\ which is given by A(c) = dim(S''^) — 
dim{Sfg^) , is also m-Lipschitz. 

Proof Here, Sioc ~ J2a ^a, where Sq, is the subset of S 
which acts trivially on a. The dimension of is given 

by 

dimiSf^^) = Rank(^H^,H5ciH<i), 

a 

where Hq, is the projector onto the set of Paulis that act 
trivially on a and H5C1 is the projector onto S"^^ . If we 
let d{ci, C2) = I we see that 



|dim(Sf„^J-dim(SfJJ| 
= I Rank(^ HqHsci H^) - Rank(^ HsHs<=2 H^) 

a a 

< Rank ^HfiHgciHa -^H^^Hse^Ho 

_ a a 

= Rank ^Ha[Hsei -Hse^lHa 

< ^Rank[(ci-C2)(Hs„)] 

a 

< m[2 Jia — dimKer(ci — C2)] = md(ci, C2). 

a 

a 

The following theorem shows that the average number 
of GHZ states which are LU-extractable from a random 
tri-partite stabilizer state is quite small, as long as none 
of the systems is larger than the other two combined. 
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Theorem IV. 5 Let \ipABc) be a uniformly distributed 
tripartite stabilizer state with local dimensions such that 
nA + fT-B > nc.riB +nc > riA, and UA + nc >nB- Then, 
the expected number of GHZ states that can be extracted 
from \ipABc) is quite small. In particular, 



ABC) 



< 



nc 



ub 



HA 



2riA+nB-i'ic 2"-4+"c^"-B 2"-B+"C-"A ' 



Proof We first express the dimension of Sioc using the 
inclusion-exclusion formula: 



dim(5^+5^+S'A) = 



dim(S'^) + dim(5'^) + dim(S'(j,) 
- dim(S'^ n S^)- dim(S'^ n S^) 
-dim(5^n5^). (5) 

Now note that the bipartite entanglement of the state 
with respect to the A\BC partition is ns + ric — dim S^, 
which must be no larger than nA- Making a similar ob- 
servation for the AB\C and AC\B partitions and adding 
the resulting inequalities gives 

2(n^-fnB+nc)— dimS*^— dim5^— dimS'p < UA+riB+nc 
so that 



dim Sj^ + dim S ^ 



d\mS^ > UA + nB + nc ■ (6) 



In order to understand the behavior of the dimensions of 
the form dim Sj^ fl 5^, first let be a fixed stabilizer on 
ABC , and Tq,Ti also be stabilizers on ABC of the same 
size as . The Clifford elements, cq, such that To = 
coS''^Cq can be put in one-to-one correspondence with the 
ci such that Ti = ciS'^'cj by using some (fixed) cq^i such 
that cq^iTqcJ^j^ =Ti, so that a uniform distribution of 
c on the Clifford group induces a uniform distribution 
over stabilizers of a fixed size for cS^c^. Thus, letting 
S — cS^^c^ with c uniform on the Clifford group, we have 

= E[#{s e S\s es^kse S^}] 



= E 



^ 1 [s e Ta & s e T£ 



.sGS 



1 [csqc^ e Ta csqc^ e Tb] 



1 [csqc^ ^Ta Sz csqc^ G Tb] 
1+ E [l [csqc^ e Ta & csoc^ e Tb] ] 

1 + (#5" - 1)E [1 [csoc^ eTAk csoc^ € Tb]] 



where we have let sq be some fixed non-trivial element 
of S*" in the last equation. 



A nontrivial Sq generates a stabilizer of dimension 1, so 
that the argument above implies that csqC^ is 
distributed uniformly on the nontrivial Paulis, leading 
us to conclude that 

E[#54 n s^] 

= 1 + (s e Ta\s e Tb,s^i) p {s e Tb\s^i) 



_ /O'ftyl+nB+'ftC 



= (2' 



■-1 



- + 1 = 1 + 



ATlA+nc ^ 



-I 



/^nA+ns+nc — \ 
1 



-I 



2nA+nB-nc ' 



where Ta (Tb) denotes the subgroup of VnA+UB tha-t is 
trivial on A [B). 

Since 4I^S ^ ^ ^b a,t least 1 and must be a multiple of 
2, this implies that P [#S^ n 5^ = l] > 1 - ^„^+i^_„^ , 
which in turn implies 

EdimS": f^Sf,< ^ . 

In a similar way we can bound E dim S'^ H S*^ and 
E dim Sj^f] Sfj to find 



E dim S^r]Sg+E dim S"^ n S'p + E dim S'^ n S'^ 



is no larger than 
nc 



riB 



UA 



which can be combined with Eq. ([n]),Eq. (O and Eq. ([3]) 
to give 



EA(S') = E(dimS'-dim5;oc) 
nc , nB 



< 



nA 



2nA+nB—nc 2"^+"C — "B 2"B+nc— "yl 



□ 



Theorem IIV.5I can be combined with Lemma IIV.4I to 
obtain the following theorem. 

Theorem IV. 6 Let \ipABc) be a uniformly distributed 
tripartite stabilizer state with local spaces of nA = an, 
nB = Pn and nc = n qubits with a, I3>1, a + l>l3 
and /3+ 1 > a. Then, the number of GHZ states that can 
be extracted from I'll'ABc) *s quite small. In particular, 
letting 5 = max(a — a-|~l) 



/ n 
p[A{ijABc) > (a + /3+l)^ + 

< exp 



64 \9{a + f3 + l) 



V. MIXED BIPARTITE STATES 

As a rule, the entanglement properties of mixed states 
can be quite difficult to understand. A mixed stabilizer 
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state, however, is always LU-equivalent to a tensor prod- 
uct of EPR pairs and a separable state, which dramati- 
cally simplifies the picture. Much like in the pure state 
case, the entanglement of a mixed stabilizer state can be 
characterized entirely in terms of the structure of its sta- 
bilizer group. The characterization we will need is given 
by the following theorem, which we will immediately use 
to get an estimate for the expected entanglement of a 
mixed stabilizer state. 

Theorem V.l (Adapted from Theorem 5 off^]) Let p he 

a mixed bipartite stabilizer state with quhits on Alice's 
system, ng on Bob's, and (nj^^ + ng — k)- dimensional sta- 
bilizer S . The entanglement properties of p can be char- 
acterized using S' , an extension of S which is the stabi- 
lizer of a purification of p to a system C. In particular, 
p is LU-equivalent to 



Eip) > 



1 



[Rank (pa ) — A + Rank (p s ) — ns ] 

(7) 

EPR pairs together with (possibly classically correlated) 
local states, where 



Proof In 3 it was shown that the number of EPR 
pairs between A and B that can be extracted by LU 
operations on a tripartite pure stabilizer state of full local 
ranks having stabilizer S is exactly ^ (dim -|-dim Sioc — 
dimS"). The state we are considering, with stabilizer S", 
may not have full local ranks (i.e., the rank of the reduced 
state on ^ or i? may be less than the dimension of that 
system) but a full rank state with the same entanglement 
can be constructed by having Alice and Bob discard any 
local pure states. The resulting state has local subgroup 
^'loc with dimension dim S[^^ > dim S[^^ — {tla + rig — 
Rank Pa — Rankp^) so that, noting that discarding the 
stabilizers of the local pure states changes the dimension 
of and S by the same amount, 

E{ps') = ^(p^,) = ^(dini5;j + dim5;'„,-dim5') 

> - (dim S'^ + dim S'l^^ — dim S' ) + 

-(RankpA + Rankps — nyi — fis)- (8) 

□ 

Theorem V.2 Let pgn-k be the a rank 2^ stabilizer state 
on AB with stabilizer 5"^'"', where A contains ua = n-\- 
alogri qubits , B has ng ~ n qubits and k ~ I3n with 
< /3 < 2 . If S is uniformly distributed. 



1 



n+-\ogn>¥.[E{ps)] > 



1 



2fe 



1 



(9) 



Proof We first consider the state whose stabi- 
lizer is 5° = {Zi, . . . , ZnA+riB-k) , where we have 
chosen some ordering of the qubits on AB. This 



state, which we call po is the reduced state on AB 
of the pure state on ABC with stabilizer S^' = 

(Z/^,Zf^^„„_,.^,ZF,Xf^^„„_,.^,XF), where i = 



UA+nB-k+j j 



1 . . . riA + ng — k and j ~ I . . . k. Letting S' — cab ® 
IcS^' c^j^g ® Ic, and Ta (Tg) denote the subgroup of 
^(nA+ufl) ttiS't is trivial on A (B), and using the ob- 
servation in the proof of Theorem IIV.5I that dim S'^ -\- 
dim S'~ -f dim 5''=, > ua + ng -\- nc as well as the fact that 
dim 5^ n 5^ = (which follows from the independence 

of {cZfBc\cZ:^X.s-k+,c\cX:^X.s-k+,cn), we find 
that 



dim5;oc > {riA + ng + fc) — dim(S'^ n S'^) — dim(5^ n S'^,) 
Considering first the expected number of elements of 

Ei^s'^ns'^ = (2"-^+"^-''-l) , +1 < 2"«-"-*-'=-^l, 



we find that 



i[dimS'^nS'^] > -logE#S'^nS'^ 

> -log(l + 2"«-"-*-'=) 
1 

> 



2n.4 + fe— "B 



and making a similar argument for S"- n S'~ we find 

B C 



Edim5;oc > riA+ng+k— 



1 



1 



Addressing the other terms in E{p) by argu- 
ing along the lines of the proof of Theorem IIII.21 
we find that ERank(pB) > ng — 2("a-"b) ^^^^ 

ERank(pA)>»a— 2"A-"B ' ^° ^^^^^ IE[£'(ps)] is no less 
than 



Ua + ng - k 



1 



1 



2nA+k+l — nB 2^^B+k+l — nA 

2 



riB - riA 



k 1 
> n TT-r 

2 2'^+! 



1 



The upper bound is obtained by noting that in 
the expression for the EPR rate given in Eq. ([5|), 
dim 5^— dim 5"= — 2k while nA+ng-\-k>6!mi S[^^. □ 

Now that we have an estimate for the expected en- 
tanglement, we would like to understand the deviations 
from this expected value. Once again, since the entan- 
glement of a stabilizer state is a smooth function on the 
Clifford group, bounds on these deviations are essentially 
immediate. 
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Lemma V.3 As a function of c € C„^+„g, the lower 
bound in Eq. ^ for the entanglement of a rank 2*^ stabi- 
lizer state PcS^ct ^-Lipschitz with respect to the metric 
defined in Eq (QJ). 

Proof That dim 5;^^ is 3-Lipschitz is immediate from 
Lemma IIV.4I Since Rank pA is simply the bipartite en- 
tanglement with respect to the A\BC partition, and sim- 
ilarly for Rank ps , they are each 1-Lipschitz by Lemma 
lin.3l As a result, the entire right hand side of Eq. ([7]) is 
■ = |-Lipschitz. □ 



3 + 1 
2^2 



which immediately implies, using the requirement 

I5iTr-e(l-"/3/2) ' ^^^^ 



P[Eips)>il + e)in--) 



< exp 



ne\l-f3/2y 



2n 



25-128 2n + a\ogn 
Combining Eq. (fTTj) and Eq. pO|) completes the proof 



11) 



□ 



Theorem V.4 Let pgn-k be the a rank 2*^ stabilizer state 
on AB with stabilizer S"^^^ , where A contains ua = 
n + alogn qubits and B has n qubits, and k = (3n with 
< /? < 2. If S is uniformly distributed, then for n suf- 
ficiently large that n > max( ^^^_^^y^^ , ^ log( ^^.^J^^y^^ )) 

and 1-^2— > max(^-^ °'ai^\ ), we have 

logn — ^ 2p ' e{l — p/2) 



P ( Eips) i il±e){n--) 



< 2 exp 



^(1-/3/2)2 



2n 



25-128 2ri,-HQ;logn 



Here E{ps)^ (l±e)(n— |) is a short-hand for the union 
of {Eips) < (l-e)(n-|)} and {Eips) > (l+e)(n-|)}. 



Proof of Theorem |V31 Lemma |V3] and TheoremUTT] 
immediately imply that 



P{E{ps) <EEips)-S) <exp 



64(riA + nB)25 



where the fact that E{ps) is |-Lipschitz rather than 1- 
Lipschitz leads to the extra factor of (|)^. This implies, 
together with Theorem IV. 21 that 



P Eips) <n ; 71-(1 - -) 



< exp 



n 

2 2* 
ne\l~P/2)^ 



2n 



25-128 2n + Q;logri 
which, using the conditions on n, gives us 



P [ Eips) <il-e) 

ne2(l-/3/2)2 



n 

2 



< exp 



2n 



25-128 2n + Q;logn 
Considering deviations above the mean, we find 

P (Eips) > n - I + I logn + n|(l - f3/2) 



(10) 



< 



exp 



ne2(l-/3/2)2 



2n 



25-128 2ri + alogn 



VI. DISCUSSION 

We have presented a general method for bounding 
the deviations of random variables which are naturally 
defined on the Clifford group. As an illustration of this 
method, we characterized the typical entanglement in 
several sorts of random stabilizer states. We found that 
a random pure state has entanglement within a fraction 
(1 — e) of the maximum possible value with probability 
exponentially close to 1 in the number of qubits being 
considered. Similarly, a random mixed stabilizer state 
with rank 2'^ and local dimensions roughly 2" has 
entanglement which is within a factor of (1 ± e) of n — -j. 
This is maximal in some sense, since it is exactly the 
entanglement of a state with the same local dimensions 
and rank which is the tensor product of a maximally 
mixed state and EPR pairs. Finally, we showed that the 
average number of GHZ states that can be extracted via 
local unitaries from a random pure multipartite state is 
close to zero and that significant deviations from this 
mean occur only with exponentially small probability. 

These results raise several questions about high 
dimensional states. Comparing the typical entangle- 
ment of a stabilizer state with that of a mixed state 
distributed according to the unitarily invariant measure 
considered in [l^ reveals qualitatively different behavior. 
In particular, a typical stabilizer state of rank 2*"' and 
local dimensions roughly 2" will have Ep = Ed = n — ^, 
whereas a typical unitarily invariant state with the 
same rank and dimensions will have Ep « n and 
Ed < — f and perhaps Ed ~ 0. It is essential to 
understand what gives rise to this difference, since it is 
the highly entangled nature of the subspace associated 
with such a randomly chosen state which makes several 
communication protocols possible [1, [iB, [3- It would 
also be nice to know whether states generated by the 
random circuit model of [l^ . [l^ have typical behavior 
more like random stabilizer states or Haar-distributed 
states. It seems clear that the behavior will be more 
like Haar distributed states when the number of gates 
in the circuit is allowed to grow exponentially [l^, but 
it would be interesting to know what happens for more 
moderately sized circuits (see also (23j). 
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The dearth of stabihzer states from which a significant 
number of GHZ states can be LU-extracted seems to be 
related to the fact that random stabihzer codes often 
fail to achieve the capacity of a very noisy channel 
[Til , m, [2^. The point is that there exist stabilizer 
codes which allow encoded Bell pairs to be transmitted 
with fidelity close to 1 in a very noisy regime where 
the average fidelity achieved by a random stabilizer 
code is bounded away from 1. These "non-random" 
codes contain states which are LU-equivalent to a large 
number of GHZ states, which explains why the codes 
are in some sense atypical. In this case, they also 
have the atypical property of allowing transmission for a 
range of noise parameter in which a typical code does not. 

Finally, we believe Theorem IH. 1 1 could be quite useful 
in the analysis of stabilizer codes. At the very least, due 
to its generality, Theorem lH.il allows one understand the 
typical behavior of a large class of random variables on 
the Clifford group without resorting to the often nasty 
computations of higher moments that would otherwise 
be necessary - given the expectation value, one need 



only compute the function's Lipschitz constant (which is 
typically quite easy) to immediately get an exponential 
bound on the probability of deviations from the mean. 
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